Thermo-statistical description of gas mixtures from space partitions 
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The new mathematical framework based on the free energy of pure classical fluids presented in 
[R. D. Rohrmann, Physica A 347, 221 (2005)] is extended to multi-component systems to determine 
thermodynamic and structural properties of chemically complex fluids. Presently, the theory focuses 
on D-dimensional mixtures in the low-density limit (packing factor r\ < 0.01). The formalism com- 
bines the free-energy minimization technique with space partitions that assign an available volume 
v to each particle, v is related to the closeness of the nearest neighbor and provides an useful tool to 
evaluate the perturbations experimented by particles in a fluid. The theory shows a close relation- 
ship between statistical geometry and statistical mechanics. New, unconventional thermodynamic 
variables and mathematical identities are derived as a result of the space division. Thermodynamic 
potentials fin , conjugate variable of the populations Nu of particles class i with the nearest neighbors 
of class I are defined and their relationships with the usual chemical potentials fit are established. 
Systems of hard spheres are treated as illustrative examples and their thermodynamics functions are 
derived analytically. The low-density expressions obtained agree nicely with those of scaled-particle 
theory and Percus-Yevick approximation. Several pair distribution functions are introduced and 
evaluated. Analytical expressions are also presented for hard spheres with attractive forces due to 
Kac-tails and square-well potentials. Finally, we derive general chemical equilibrium conditions. 

PACS numbers: 02.50.-r, 05.20.Jj, 34.20Cf, 51.30.+i 



I. INTRODUCTION 

The detailed knowledge of atomic and molecular po- 
pulations is essential to analyze and evaluate accurately 
thermodynamic properties and chemical processes in flu- 
ids as well as their monochromatic opacities. In partic- 
ular, they are required in studies of stellar interiors, at- 
mospheres and circumstellar envelopes. Gas densities in 
stellar atmospheres and circumstellar structures are low, 
they are hardly higher than 0.001 gem -3 . Nevertheless, 
in these cases it is necessary to deal with a high variety 
of particles and electronic configurations and therefore, 
with a series of physical processes to calculate the radia- 
tion emitted by the star accurately Q, . 

Among the most difficult challenges in determining gas 
models is to obtain self-consistently atomic densities and 
opacity absorption coefficients by taking into account the 
non-ideal effects. The source of non-ideal effects is the 
presence of Coulombic interactions between charged par- 
ticles as well as due to the short-ranged interactions be- 
tween neutral-neutral and neutral-charged particles. Due 
to its relevance and complexity, the elaboration of gas 
models from mechanical statistical approaches is still an 
active area of research today [3,OLIHISl1IIIE1E3>EI>0]- 

The present study is motivated by the difficulties of 
predicting quantitatively the optical properties of stel- 
lar atmospheres. In typical gas models the opacity data 
are obtained through a complimentary calculation once 
the equations of state and mean abundances of bound 
states are obtained. This procedure can produce, how- 
ever, internal inconsistencies in a model, mainly when 
the abundances of atoms are calculated using internal 



energies that are different from those assumed to obtain 
the radiative transitions 0, . This can lead to seri- 
ous disagreements between predictions and observations 
or experiments (e.g., the unphysical Lyman opacity |15l|). 
Therefore, it is highly desirable to develop a theoretical 
formalism that treats atoms subjected to different per- 
turbations by considering them in separate groups and 
evaluate accordingly the populations and opacities in a 
consistent way. Such a gas model has then to consider 
explicitly not only all atomic excited states but also do 
they need to discriminate different perturbation states 
with the use of an appropriate particle-state variable. 
Our work is then oriented to develop a model with these 
characteristics. 

With the aim of establishing a detailed and accurate 
equation of state for fluids, in a previous work we have 
developed a gas statistical formalism which combines 
free energy minimization methods with space partitions 
[l3|. The novel feature of this formalism, as compared 
to other more standard and well-established statistical- 
mechanical theories of the fluid/liquid state, is the at- 
tempt of keeping trace, in a thermodynamically self- 
consistent way, of a more refined structural information 
on the local environment of a given particle through an 
extra parameter, the "available volume" , which is deter- 
mined by the distance of the reference particle to its clos- 
est neighbor. For noninteracting particles, the available 
volume 'V is equivalent to the spherical volume whose 
radius is the distance between the respective centers of 
the particle and its nearest neighbor. This volume is 
considerably smaller for particles with repulsive interac- 
tions. The theory [l^ | was formulated for one-component 
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dilute fluids, where the Helmholtz free energy is writ- 
ten in terms of the occupation number distribution N v 
of the w-variable and where the particle interactions are 
introduced using pair potentials. The result thus ob- 
tained represents a unified treatment of thermodynamics 
and the structure of fluids. The formulation was applied 
to hydrogen atoms in an electrically neutral medium of 
charged particles encompassing ions and electrons. The 
atomic populations were derived by establishing groups 
of atoms having different plasma perturbations measured 
according to the size of their available space volumes v. 
Although the results obtained are only illustrative of its 
possible performances, this theory represents the first at- 
tempt of evaluating in detail the populations by combin- 
ing consistently the plasma interaction effects. 

Approaches based on space partitions were attempted 
in the past with the so-called "cell" theories to obtain 
equations of state for fluids. However, the character of 
these space divisions is quite different as compared to 
that em ploy ed in [131 and in the present paper. The cell 
models |l6l \vH Il8[ were originally developed by Eyring 
and Hirschfclder pf9| and by Lennard-Jones and Devon- 
shire [2(J. They were inspired by the lattice-like config- 
urations acquired by some systems at high densities. In 
such theories the fluid volume is divided into an imagi- 
nary lattice of cells. Each particle is then confined into 
a cell, which implies that a space partition is assumed 
beforehand, instead of being deduced from the model 
proper. Besides, by ignoring the interchange of parti- 
cles between cells, an error is introduced in the system- 
entropy calculation, which is usually corrected using an 
arbitrary "communal entropy" factor. 

The perspective of the theory presented in |13| is new, 
as it determines a space division by assigning the v- 
volume to the particles. This is achieved by minimizing 
the system free energy under the condition that the sum 
of all v- volumes be equivalent to the entire volume of the 
fluid. In this procedure, no particle confinement is estab- 
lished, because the assignment of an individual volume 
to each particle does not introduce any ownership con- 
nection between space regions and particles, i.e., v does 
not carry any information on the particle coordinates. 
Therefore, the formulation presented in [l3j is free from 
difficulties that are inherent to the cell theories, i.e. en- 
tropy corrections, and it can provide new results, such 
as the combined thermodynamical and structural infor- 
mation of a fluid, which cannot be derived from the cell 
models. 

In the present paper we generalize the formalism to 
D-dimensional multi-component fluids that was initially 
developed in [l3| only for pure or single-component sub- 
stances. It aims at giving particularly simple analytic 
forms to provide a unified derivation of thermodynamic 
and geometric properties that characterize them. The 
formulation of the theory is general, but attention is fo- 
cused on its application to hard sphere models, which 
represent the repulsive interactions among atoms and ge- 
ometrically simple molecules [2l]. We also include the 



analysis of hard spheres with weak long-range attraction 
(Kac type) and hard spheres with square-well potential. 
The approximations discussed here should be suitable for 
more realistic short-ranged potentials, like those govern- 
ing many atomic and colloidal fluids. Our introduction 
of space partitions in statistical thermodynamics marks a 
considerable progress in the description of gases with re- 
spect to previous theories, since it gives the possibility of 
obtaining thermodynamic and structural information on 
fluids simultaneously and self-consistently. The present 
extension to multicomponcnt fluids enables us to treat 
arbitrarily complex chemical mixtures by including sys- 
tematically the inter-particle interaction potentials. The 
theory can have many applications, not only in Astro- 
physics, but also in other fields such as the analysis of 
colloidal and protein solutions in biochemistry. 

In the following sections, we first introduce a descrip- 
tion of states of a gas mixture based on space partitions. 
General expressions for the Helmholtz free energy are 
presented in Sect. 1 1 1 1 1 and the thermodynamic equilib- 
rium states are derived in Sect. IIVI We analyze there 
some non-conventional thermodynamic variables that are 
generated by the theory and derive general expressions 
to calculate chemical potentials. Then, applications of 
the formalism to a few multi-component systems based 
on hard particles with and without attractive forces are 
given in Sect. [V] The results obtained are compared with 
those from Percus-Yevick and scale particle theories. In 
Sect. IVII we extend the theory to include chemical reac- 
tions in the gas. The conclusions are presented in Sect. 

eh 

II. GAS STATES 

The approach developed here is inspired by the vari- 
ety of spatial patterns shown by fluids where the particles 
have different types of interactions. The main assump- 
tion that will be made is that the spatial structure of a 
gas in equilibrium can be described by an appropriate 
division of the space among the particles. The study of 
single-component fluids at low densities ^| has shown 
that the volume v assigned to a particle depends on the 
interactions the particle has with the closest neighbor. In 
Sect. El we give a formal definition of the closest neigh- 
bor. According to the mentioned result, the generaliza- 
tion of the formalism to multi-component fluid systems 
(mixtures), requires that the particle population Ni of 
each chemical species i be divided into sub-groups or 
sub-classes in accordance with the chemical species of 
the nearest neighbor, i.e. 

N i = N il + N i2 + N i3 + ..., (1) 

where the sub-indices 1,2,3... specify the various chem- 
ical species of the closest neighbors. We consider that 
for any reference particle in whatever configuration of N 
particles, there is 'one and only one' jth neighbor, with 
1 < j < N — 1. This means that if there is degeneracy 
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where two particles can be equidistant from the refer- 
ence, one of them will be designated as the jth neighbor 
arbitrarily, while the other will be the (j + l)th neigh- 
bor. In what follows we shall use the sub-index ilv for a 
particle of class i that has individual available volume v, 
whose nearest neighbor is of class /. The state of the gas 
will then be described by the following collections: 



N nv N 12v 
N 2 iv N 22v 



(0<v< V), 



(2) 



where Nu v dv is the number of particles ilv with avail- 
able individual volumes ranging from v to v + dv. The 
physically acceptable states in J3J must obey the particle 
number Nu and gas volume V conservation conditions 
given by 



Na 



^ = E / vN avdv. 



(3) 



(4) 



Thus, the individual volume v assigned to each particle 
acts as an attribute inherent to each particle and the 
sum of f-spaces is equivalent to the entire volume of the 
gas. This space partition does not involve any knowledge 
of particle positions, since the assignment of a v-volumc 
to each particle does not carry any information on the 
particle coordinates. 



have then that the total number of ways W to generate 
a set {Nn v } of occupation numbers from the individual 
complexions W% is W — J\ i Wi, where: 



(dv/V) N «° dv 
(N ilv dv)\ 



(8) 



Therefore, according to [l3( the configuration entropy is: 



S conf = klnW = -kJ2 I 
si Jo 



V N av hx(^-)dr. (!)! 



which finally determines J7J), since F con f = —TS con f. 

The interaction can be generalized in the same way as 
follows 



U,, 



E 

a 



2V 



dv, (10) 



where the pair-interaction factors 



are: 



V' 



uav,j(u)gu v ,j(uj)d(j. 



(11) 



In Hip . UjXvj is the interaction potential between a par- 
ticle of class ilv and another of class j. The distribution 
function of particles j with respect to the reference par- 
ticle ilv is written as gu v ,j- We note that the integration 
variable u> in Eq. (jl 1|1 is a spherical volume centered at 
the particle of class i and that V denotes the minimum 
lu— volume able to enclose the entire volume of the mix- 
ture. The total free energy of the gas-mixture is then: 



III. HELMHOLTZ FREE ENERGY 

To evaluate the free energy of the gas we use a com- 
position of translation-, configuration- and interaction- 
dependent terms: 



F, 



conf 



Uint , 



(5) 



which is suitable for a large variety of fluids 22]. The 
first two contributions to F can be written respectively 
as: 



F h 



E 



— ln {- 



F 



conf 



E 



P 



In 



Ng v V 
Ni 



dv, 



dv, 



(6) 



(7) 



where (3 = (kT) -1 , k is the Boltzmann constant, T is the 
temperature and is the thermal wavelength of par- 
ticles of class i. F con f was obtained following the steps 
described in [l3| , except that in the present case was done 
for a gas-mixture. For the multi-component system we 



^E 



(3 U \ N, J ^ 2V 

3 



dv. (12) 



IV. EQUILIBRIUM STATES 

The equilibrium distributions Nu v are obtained by 
minimizing the free energy l|12f) and taking into account 
the conservation conditions given by Eqs. 10 and 10} . 
The minimizing procedure must take into account that 
the variations 5Nn v affect the total particle populations 
of each species according to the relation: 



N 3 = E / N -imvdv. 



We thus obtain for the state of equilibrium 

Ni 

Nuv = — £ exp [-(3 (-yv - na + e Uv + e 
\ 

where 

Nj (f>Uv j 



E 



2V 



(13) 



(14) 



(15) 
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and 



potentials fjtn, 



V — 

^ 2V 



(16) 



In Eq. 1)14(1 , [in and 7 are the Lagrange multipliers corre- 
sponding to conditions © and 10} respectively. Accord- 
ingly, the potentials /in determine the total population 
of particles i and their redistribution into sub-sets as a 
function of the nearest neighbor class. The variable 7 de- 
termines the distribution of individual available volumes 
v. 

It is useful to note that the total interaction energy 
can also be written as 



E 



E 

a 



Nu v eu v dv, 



(17) 



which enable us to see that q is the average energy due 
to all interactions the particle of class i has with the 
gas, while eu v represents the interaction energy of the ilv 
particle with the gas. Finally, using l|14l) and l|12l) we can 
write the equilibrium free energy for multi-component 
gases in the Eulerian form 



F = - 1 V + J2tMiN a -U l , 



(18) 



= E ^ h 

I III 



dN, 



irn 



dNi 



(22) 



T,V,Ni 



In fluids with low enough density, it is expected that 
the splitting up of the Ni populations into sub-groups 
Nu [Eq. ||TJ] be determined by the mole fractions xi = 
Ni/N of species present in the gas. This means that the 
population of particles i having as closest neighbors the 
class I particles is proportional to the total population of 
% — particles and to the abundance of I— class particles: 



N a = 



NN... 



N 



= xiNi 



(23) 



We can also say that the probability of a particle had 
an /—class as the closest neighbor is given simply by the 
abundance of I species. In this case (1221) becomes 



Mi = yiifMl + m) x l ~ E VlmXlXm- 
I 



(24) 



It is easy to show that if fin is given by the sum of two 
terms fi + gi, i.e., one term that depends on the class i 
and the other on I, it will be fM = fi + <?i- This is what 
happens in a mixture of ideal gases, as it is shown in 
what follows. Let us note that introducing l(Tl|l into © 
leads to the useful relation 



A. Lagrange parameters 

Relation (|A5|) obtained in (Sect. [XJ and ijTSjl help 
us to show easily that fin behaves as a thermodynamic 
potential associated to the Nu population: 



fJ-u 



dF 



dNi, 



(19) 



T,V,N jr , 



Similarly, making use of (|18fl and l|A10|) we find that 7 
relates to the pressure P in the already known way for 
single-component fluids jLj], 



P = 7 - 



Un 



V 



(20) 



i.e., the gas pressure is the sum of 7 and the interaction- 
energy density. 



B. Chemical potentials 

The chemical potential corresponding to i— class par- 
ticles is: 



dF 

m 



(21) 



T,V,N j¥ ,i 



The substitution of (|18|) into l|21|l and the use of relation 
(|A13|) lead to the following expression of fii in terms of 



Nu 



\7 D N i Z il e^ a 



where 



exp [-(3 (7« + eu v + e,)] dv. 



(25) 



(26) 



On one hand, for non-interacting gases (hereafter identi- 
fied with the superscript *), we have 



(07)- 



(27) 



which is obtained from (0} and the estimation is done for 
the thermodynamic limit (N, V) —> 00 with n = N/V = 
constant. Therefore, for a mixing of gases it becomes 



ml 
v 



(28) 



which through Eq. I|24|l yields the following exact result 

if 



Mi 



fcTln 



V 



(29) 



On the other hand, for interacting fluids, relation 1261) 
gives 



with 



(fii) 



e -(x+Peu*) dx 



(30) 



(31) 



5 



where x = f3jv. In general, if interactions in the gases 
are fairly low, factors i|31|) will be qu fa 1. In fact, qa — > 1 
if — *0. Using the condition (|23[1 one obtains 



fcT 



= In 



07 



fcT 



- In qu , 



in which 



\n(quqii) - ^ x m ln^. 



(32) 



(33) 



(34) 



The excess chemical potentials /tiff = /Lt^ — and 
fi^ — /i* are formed by three contributions with spe- 
cific characteristics. The first right-hand term of equa- 
tions i|32*|) and represents the deviation of the divi- 
sion of space among particles with respect to that in a 
perfect gas, i.e. /3j ^ n. This effect is the same for all 
chemical potentials. The second right-hand term in H32|) 
and Q33JI gives the contribution due to the potential en- 
ergy acquired on average by an z— class particle due to 
its interaction with all the remaining particles in the gas. 
The third term corresponds to interactions that affect the 
distribution of individual spaces v, as can be seen from 
definition of qu in (|31|l . This last contribution can be 
significant even for interactions that do not carry energy 
to the gas, as it is for infinitely strong repulsive forces 
(see Sect. 0. 

By virtue of Eq. I|33|l . we can obtain a similar expres- 
sion to H25|) for the total population of a given class 



Ni = A. 



- D NZie^ 



(35) 



of non- interacting particles, see Sect. [I| reduced by an 
amount a* ; that depends on the repulsion between the 
particle and its neighbor with the shortest contact dis- 
tance (hereafter nscd): 

v = v*-a* l . (38) 

Other interpretations than (|38|l are also possible. How- 
ever, the results that are obtained in the present section 
strongly support the choice given in (|38[l. As a conse- 
quence of this interpretation, a particle of class j =/= I 
with ay < an can have its center at a shorter distance 
from the center of particle i than the particle I. 

The identification of a volume per particle enable us 
to specify the form of the distribution function of pairs 
(hereafter pdf) gu v ,j{w) for low densities. Notice that 
9iiv,j(^)duj gives the probability of finding a particle of 
class j located between the surfaces of spheres to and 
to + duj both centered in a particle ilv. According to the 
results obtained in for single-component fluids, we 
deduce that the pdf with j = I has a contribution due to 
the nscd given by a Dirac delta function S (x) evaluated 
on the nscd center and divided by the density of I— class 
particles. In fluids of low density, the contribution to the 
pdf due to the remaining particles j can be approximated 
with a Heaviside step function 0(x) [Q(x) = 1 for x > 
and Q(x) = when x < 0], which implies that the 
position of these particles is uncorrelated with that of 
the reference ilv one. Thus, 



guv,j (oj) =—5(u — v — ali)+Q(u—v—a^i+au—ai 



)■ (39) 



where 



The terms in the argument of the step function warrant 
that the nscd is a I— class particle. The Kronecker delta 
symbol Sji ensures that the contribution due to nscd is 
taken into account only if j = /. Introducing Eqs. I|37(l 
and J2HJ) into Eq. ifTTl) we obtain 



Zi = G97) _1 0i 



(36) 



oo; v < b a 
0; v > b a 



(40) 



HARD PARTICLES 



and 



A mixture of hard spheres in D dimensions is charac- 
terized by interaction potentials described by: 



where 



(41) 



oo; lu < a 
0; lu > a 



2 D o~ij, where an is the volume 

di)/2, 



(37) 

of a D- 
being di 



with ftjj — ^ u ij- " "i. i 
dimensional sphere with diameter (di 
and dj the respective diameters of particles of classes i 
and j. The interaction does not depend neither on the 
I— class neighbor nor on the available volume. To iden- 
tify the available volume v of a give n particle we adopt as 
reference the results found in |13| for single-component 
fluids. We adopt then as ansatz that for low particle den- 
sities the available volume is the ideal volume v* (i.e., v 



ba — an — a*[. 



It follows then that 



oo; v < b u 
0; v > b a 



(42) 



(43) 



while ti = Vi. In the thermodynamic limit, Eqs. © 
and |g} yield 



fia = kT 



b u n t ^ / n a Xf n 



1 — bn 



nj 1 — bn 



(44) 
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FIG. 1: (Color online). The compressibility Z factor calcu- 
lated with Eq. 146H (dashed lines) and Percus-Yevick results 
pfij l (dotted lines) as a function of the packing fraction r\ for 
3D binary mixtures with x\ — 0.5 and (fe/di = 0.01, 0.4, and 
1.0 (from bottom to top). 



FIG. 2: (Color online). Excess chemical potentials of a 
diluted solution (here species 2) calculated with Eq. (|50P 
(dashed lines) and those from the scale particle theory [2q | 
(dotted lines) as a function of the packing fraction r\ for 3D 
binary mixtures with X2 = 10~ 5 and da/di = 0.5, 1.0, 2.0, 
3.0, and 5.0 (from bottom to top). 



nav = &(v - bu)na0jexp [~Pj(v - b a )] , (45) 



1-bn' 



where we have written 



bn = b. 



(46) 



(47) 



which of both, the reference or the nscd particle, is the 
largest. Nevetheless, as it is expected, a big nscd reduces 
more the space v of a given particle than a small nscd. 

On the other hand, from Eqs. (|24H and 14411 we obtain 
the excess chemical potential for hard spheres of species 
i that is given by 



kT 



ln(l - bn) + 



km - bn 



1 — bn 

which in particular for D = 1 reduces to 



(50) 



As it is expected, the gas does not store interaction en- 
ergy, which implies that [Zj„t — and thus 7 = P. From 
the comparison of Eq. (|46Jl with the corresponding in the 
Virial expansion for mixtures of hard spheres [16| , enable 
us to make the following identifications 



ba = 2 



D-l, 



era, 



nu = xini 



(48) 



(49) 



Relation (|48|l is nothing but an extension of previous 
results found for single-component fluids, while 149(1 con- 
firms the relation already established in Eq. (|23|l . Since 
= &u the quantities an, bu, and a* t are symmetric in 



the indices il. In particular, a* 



vD-l 



0,l 



so that 



the volume of a hard sphere is reduced with respect to 
the ideal value v* by an amount that does not depend on 



kT 



m l 



vim 



G vi n 



1 - Y^i °u n i 



(51) 



Eq. gSJ is exact for D = 1 jU UM an d ^ S ives om y 
an approximate solution for D > 1. However, it gives 
the correct asymptotic behavior for n —> 0. This is 
shown explicitly in Fig. ^ where the compressibility fac- 
tor Z = P/(nkT) is calculated with Eq. I|46() and com- 
pared with the Percus-Yevick (PY) approximation |2{| 
for a 3-dimensional binary mixture, as a function of the 
packing fraction n = n^au which gives the fraction of 
the total volume occupied by the hard particles. Since 
we compare our results with those of the PY theory, one 
of the most successful theoretical approaches in studies 
of fluids, let us recall that it is based on an approximate 
integral equation of the radial distribution function from 
which the thermodynamics of hard-body systems follows 
as a by product. In few cases, it can be solved analyt- 
ically, e.g. mixture of hard spheres, and the solutions 
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FIG. 3: Excess chemical potentials obtained with Eq. I50H 
(symbols) and Percus-Yevick theory results 27j (lines) as a 
function of the diameter ratio dijd\ for 3D binary mixtures 
with n = 0.02 and x 2 = 0.05, 0.1, and 0.5. 



obtained are very close to those issued from computer 
simulation [22|. 

Likewise, Eq. ifBTJl is exact for hard rods jl]} ■ Eq. fBUjl 
provides the correct low-density behavior of the excess 
chemical potential for D > 1. This last is illustrated in 
Figs. |3 and which show the comparisons of our results 
obtained for D — 3 with those issued from the scale par- 
ticle theory j2|j , which is correct for a solute infinitely di- 
luted and with the PY solutions found by Lebowitz [25| , 
respectively. It is important to note that all expressions 
derived for mixtures of hard-spheres are exact up to the 
second term in the virial expansion. Finally, we can con- 
clude that the achievement of agreements shown in Figs. 
^ 13 and |31 confirm the appropriateness of the closure rela- 
tions 1|23[) and l|38fl adopted, where a* t = 2 D ~ 1 au. Since 
our work aims at giving a detailed representation of com- 
plex gases at low densities (77 < 0.01), the results we have 
obtained can be considered as very satisfactory. 



A. Pair-distribution functions 

From detailed ga v ,j functions it is possible to calcu- 
lated the averaged pdf with the occupational numbers 
nu v (see 0). For example, the probability of finding a 
j particle between the spherical surfaces oj and lu + dto 
centered in a i particle, regardless of its available vol- 
ume, whose nscd is an I particle, defines a conditional 
pair distribution function cpdf which is given by 

1 f°° 

9Uj(w) = — / gu v j(uj)nu v dv. (52) 

n il JO 



The traditional can be computed from the summation 
over all species of cpdf appropriately weighted by the 
probability of finding a nlcd of a given class, 



9ij 



H =y ^xiguj{u). 



(53) 



For mixtures of hard-spheres, Eqs. 1521) and (|53|) give 
respectively, 



gu,j(u) = Q(u> - aij) 



U i; ^-i: 



g ij (w)=®(u-a ij )n + - 



bn 



bn 



exp 



1-bn 



(54) 



(55) 



Eqs. I|54l) and (|55() are approximations suitable for low 
densities, because the function gu v ,j in Eq. Ij39j) takes 
into account explicitly only the correlations between the 
reference particle and its nscd. However, they contain the 
most important deviations with respect to the solutions 
for perfect gases. 

To our knowledge, this is the first time that the cpdf 
are suggested and calculated for mixtures of hard spheres. 
Figure 01 shows some cpdf for a binary mixture with 
i] = 0.1 and X\ = 0.7. In the chosen example, particles 
of class 1 have a diameter twice as large as those of class 
2. Functions gu.j become zero for distances r between 
particles i and j shorter than the contact distances, as it 
happens for the conventional pdf g^. Nevertheless, con- 
trary to functions gij, the behavior of cpdf at moderate 
distances can be strongly affected due to the informa- 
tion on the nscd class carried by the reference particle. 
Thus, up to the present level of approximations, the func- 
tion <7i2.2 grows above 1 for small values of r (r> 1.5i?i, 
r— ►1.5i2i, where R\ is the radius of class 1 spheres). This 
is due to the fact that the probability of finding a particle 
of class 2 at a distance r w 1.5i?i from another of class 1 
augments if it happens that this last has a nscd of class 2. 
This augmentation concerns the non-conditioned proba- 
bility of finding a class 2 particle and so, it is higher the 
lower is the abundance of the considered species. Con- 
versely, for the same range of distances, the function 311,2 
decreases because the possibility of nscd be of class 2 is 
removed. 

The results corresponding to a mixture of ideal gases 
are contained naturally in the expressions l|54|l and 15511 . 
In this case, the results are exact. For permeable spheres, 
i.e. = bij = a*j =0, the expected correct result g*j = 1 
V i, j is obtained. Moreover, 



1 - e~ nu 
1 + {x- 1 - l)e 



(7=0 



(56) 



Figure [S] shows some cpdf corresponding to a binary bi- 
dimensional mixture predicted by Eq. (|56|l (full lines), 
which arc in excellent agreement with the results ob- 
tained from numerical simulations (dotted lines) based 
on two sets of uncorrelated points created by a random 
number generator. 
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B. van der Waals mixture 
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FIG. 4: (Color online). Conditional p.d.f. from Eq. JMJ as 
a function of the normalized radial distance r/Ri (Ri = 4p is 
the radius of a particle of species 1) for a 3D binary mixture 
with rj = 0.1, x\ = 0.7, and dij&\ = 0.5. 
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FIG. 5: (Color online). Conditional p.d.f. computed with 
Eq. (full lines) as a function of the dimensionless ra- 

dial distance r^frvi (in units of the mean distance between 
particles of species 2) for a 2D binary mixture of perfectly 
permeable particles at x\ — 0.7. The results obtained from 
numerical simulations with 10 Poissonian points are shown 
for comparison (dotted lines). 



The van der Waals model for a mixture can be obtained 
by adding to the repulsive forces [Eq. l(H7jl] a generalized 
Kac type attraction potential, 



Uil v j(u)) 



00: uj < a 



(57) 



where > are constants and J is a function that 
satisfy the condition J °° J(uj)duj = 1. In the limit 7*^0 
it becomes, 



'ilvd 



00; v < bu 
-en; v > b a 



and 



so that 



- 2 J2j e '. 



00; 



v < bu 
v > b a 



(58) 



(59) 



(60) 



The attracting forces modify the chemical potentials and 
the total interaction energy of the gas, so that it results, 



fJ.ii = kT 
and 

where 
and 



nb, 



1-nb 



-In 



n,X D 



rii 1 — nb 



Hnt 



(61) 



(62) 



(63) 



(64) 



Nevertheless, Eqs. (|45H and (|46|l remain unchanged as 
it happens for single-component fluids [l3|. We then see 
that using Eq. l(2Tj|) the known relation P — (nkT)/(l 



nb) 



2 /2 is recovered. Relations (|63f) and (|64l) are 



often used as mixing rules which relate the properties 
of the pure component models to those of mixtures in 
extended models. They are considered adequate for mix- 
tures of monopolar compounds p9T l30l | and vapor-liquid 
phenomena [3l|. 



C. Attractive hard spheres 

A simple model of hard spheres having short-range 
attractive forces, as opposite to the long-range van der 
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Waals model, is obtained by adding a square-well poten- 
tial at the edge of a repulsive nucleus, 



{oo; lo < a 

—6ij; a tj < u) < aij + £ , (65) 
0; ijj > + £ 

where £ is a characteristic attractive volume surround- 
ing the ay volume of the attractive nucleus. For sim- 
plicity we adopt a value of £ independent of the con- 
sidered pair of particles. This potential is used often to 
model the interactions among colloidal particles and pro- 
tein molecules j33, HH, Q , which commonly have attrac- 
tive forces in the very nearest vicinity of their surfaces 
(C<Oij). With Eq. (jnni, Eqs. and fiB lead to 

oo; v < b a 

'■ir.i = \ _e w(tf + b u+£~v) 5 ki<v<ba+£ , (66) 

0; v>b a +Z 



and 



oo; v < b 



where 



-tf-Bi(bu+t-v); ba<v<bu+£ , 



0; v > bit + £ 



6ij rij 



(67) 
(68) 



The state equations that result from the particle and vol- 
ume conservation relations and are: 



Hu = kT 



p-ybu + hi 



Xu 



where we have 



-Pit . 



0(e u /2+B^) 



-P(~/+b)t 



Pil + Bi) 



(69) 
(70) 

(71) 



Y a = [1 + /3 7 {hi + 0] 

+ (A) 2 e^ 2+B ^ {1 + /3( 7 + (72) 
- [1 + /3( 7 + Bi)(6« + 0] e-^(T+ B «)C} . 

Besides, the ^-occupation number distribution is given 
by 



Xu 



(73) 



We can also obtain an analytic expression for the inter- 
action energy: 



Uint 
V 



V- Pfnu 0( ei ,/2+B i £) 



+B 



l-[l+/3( 7 +a,k]e- g(7+B ' K 
[Ph+B, 



(74) 



)}• 



It can be interesting to have limiting expressions of these 
equations for small enough values of £. Keeping only the 
zero and first order terms, we derive 



X a = 1 



■Pit (' 



,/W2 _ 



Y, 



(75) 



(76) 



If also the potential wells are the same, en = e V(z, I), the 
state equation (|7U[) reduces to a quadratic relation in /?7, 



(/3 7 ) 2 ^ £/2 -l)+/37 



1 — bn 



(77) 



where b = J^u %i%lbu- When the thickness £ or the depth 
e of the potential well tend to zero, the state equation lf??|) 
recovers the form corresponding to a mixture of pure hard 
spheres. The same result is obtained for high enough 
temperatures, kT ^> e. 



D. Space distributions 

The thermodynamic variable 7 has the function of ad- 
justing the allotment of space among particles in order 
to satisfy the mathematical requirements of space filling 
established by In mixtures of non-interacting par- 
ticles the occupation number distributions nu v depends 
on the individual volume per particle v through the ex- 
ponential exp(— fiyv). In this case, the scale factor f3-f is 
equivalent to the total number density of particles n. Ob- 
viously, the allotment of space is not homogeneous. Most 
particles have smaller available volumes than the average 
volume per particle, n _1 , while the statistical partition 
of big volumes is low. 

The particle interactions change the assignment of the 
volume v for each particle and thus they change the oc- 
cupation numbers nu v , as compared to those for ideal 
gases. Therefore, the interactions also affect the 7 val- 
ues. In Fig. [fJlwe have plotted fi^n -1 as a function of the 
packing fraction 77 for perfectly permeable spheres (dot- 
ted line), pure hard spheres (full line), and hard spheres 
with square- well potentials (dashed lines). Notice that 
fl'yn is a measure of 7 in units of perfect gas values. 
Figs. \7\ and [S] show occupation number distributions nu v 
for binary mixtures of pure and attractive hard spheres, 
based on Eqs. and lf75)l. respectively. The lower cut- 
off value for each curve corresponds to the smallest vol- 
ume v permitted by the repulsion between a hard sphere 
and its nscd, as given by Eq. i|38[l . The square- well po- 
tential yields the peaks observed in Fig. [S] If we make ab- 
straction of the normalization constants - area under the 
curves representing the nu population - the peak height 
of curves nu v vs. v in Figs. and |H1 depend on the 
abundance of the corresponding pairs il. In the frame 
of models of mixtures considered here the distributions 
ni2 V and n-iiv are coincident. 
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FIG. 6: Factor /^yn -1 as a function of the packing fraction 
for 3D binary mixtures of spheres with xi = 0.7 and cfe/di = 
0.5. Dotted line: perfectly permeable spheres; full line: hard 
spheres and van der Waals fluid; dashed lines: square-well 
fluid with £/<t = 0.2 (a is the mean volume of a particle) and 
reduced temperatures kT/e = 0.12, 0.5. The result for pure 
hard spheres is obtained for kT/e — » oo. 



In a hard sphere mixture, the repulsions among the 
particles produce an increase of the mean volume per 
particle. As a consequence, at given temperature and 
density the value of 7 must increase (full line in Fig. [BJ 
to warrant the occupation distributions nu v be correctly 
displaced towards higher values of v (full lines in Fig. 
0) . On the contrary, the attractive forces of finite range 
in a square-well fluid, increase the population of parti- 
cles with small individual volumes (full lines in Fig. |8J) 
and reduce the values of 7 as compared to those of pure 
hard spheres (dashed lines in Fig. [fJJ). These changes de- 
pend on the temperature value. Low temperatures favor 
the particle adherence and thus increase the occupation 
number of low available volume. 

The long-range interactions of Kac type analized in 
Sect. IV Bl are such that the magnitude of the interac- 
tion of each particle with the remaining medium does 
not depend on the nearness of its neighbors. Hence, the 
attractive forces have no incidence on the space occupa- 
tion distribution nu v neither on 7, both being equivalent 
to those in a mixture of pure hard spheres and hence 
independent of the temperature. 

It is useful to point out that, for ideal and pure hard- 
sphere mixtures, the quantity (3^n~ l coincides with the 
compressibity factor Z = P/nkT since these gases cannot 
store interaction energy and, therefore, 7 is equivalent 
to the gas pressure P. It is not true for the van der 
Waals and square-well models, because they have a non 
negligible interaction energy, so that 7 and P differ as 



FIG. 7: Dimensionless double densities nu v /n 2 as a function 
of the available volume per particle (in units of volume o\\ 
of a particle of class i = 1) for a 3D binary mixture of hard 
spheres with r\ — 0.1, x\ = 0.7 and dijd\ = 0.5 (solid lines). 
Results corresponding to fully penetrable spheres are shown 
by comparison (dotted lines). 



described by relation l|20|l. 



VI. CHEMICAL EQUILIBRIUM 

In this section we extend the theory somewhat to in- 
clude cases where chemical reactions can occur. In this 
case, the total free energy contains a contribution given 
by the summation of all particle internal energies, 



E = 



.-I Jo 



NuvEidv, 



(78) 



where Ei is the internal energy of a particle i and the ref- 
erence zero energy is considered the same for all species. 
If the internal energies Ei are degenerate with multiplic- 
ity Wi, then Wt in Eq. (JSJ becomes, 



Wi = A^n 



(wjdv/V) N ^ dv 
(N llv dv)\ 



The configuration free energy takes thus the form 



F, 



conf 



E 







In 



NuvV 
WiN t 



dv. 



(79) 



(80) 



The results obtained in previous sections given by (|14l 
I5U1 can be generalized now as, 

Nilv = -rh w i ex P [-0 il v - Ml + e Uv + £ * + E i)] ) ( 81 ) 



11 



1.5 

c 

^ 1 
a 

0.5 



1 


1 


1 1 


1 1 1 1 1 1 1 

kT/e = 0.5 ~ 




(il)=(ll) 


£/cr = 0.2 1 
77 = 0.1 ~ 


(12) 




- 


1(22) 
1 






~r~, — 1- 1 





1 


1 1 


5 10 



FIG. 8: Dimensionless double densities nu v /n 2 as a function 
of the available volume per particle (in units of volume an of 
a particle of class i — 1) for a 3D binary mixture of attractive 
hard spheres with r/ = 0.1, x\ = 0.7, ab/di = 0.5, £/<r = 0.2 
and kT/e = 0.5 (solid lines). Results for pure hard spheres 
are shown by comparison (dotted lines). 



that gives the generalized form of the mass-action law. 

The results given in this section bring useful tools to 
study gaseous mixtures, including their chemical equilib- 
rium. Eq. (|81(l sums up the main feats of the present 
formulation. It enables to describe populations in react- 
ing gases at a level of details attained never before, i.e. 
to evaluate populations of atoms and molecules in a con- 
sistent thermodynamic way, not only by obtaining the 
density of each species, but giving as well its populations 
per groups according to the available volume of each par- 
ticle, where this variable is linked directly with the de- 
gree of perturbations experienced by the atoms due to 
the presence of neighbors. In the preceding work |13| we 
have foreseen some results that could be obtained with 
the present theory. For example, it is possible to show 
that atoms of a given species in a gas do not have all the 
same effective ionization potential (electronic energy rel- 
ative to the energy threshold corresponding to free elec- 
trons), but can they have a variety of values according to 
the available volume of particles. The application of the 
present theory to mixtures of atomic particles with more 
realistic interaction potentials than considered here will 
be presented in a forthcoming paper. 



VII. CONCLUSIONS 



and 



(82) 
(83) 



Now, let us consider the reaction 

aA + bB + ...< — >cC+..., (84) 

that involves the species i = A, B, C, ... and the stoichio- 
metric factors a, 6, c, .... The equilibrium condition of the 
reaction is 

dF = zZ\^W )dNi = J2iHdNi = 0, (85) 

. \".\, TA.X . J 



which is subject to the conditions 

dN A dN B dN c 



(86) 



a b c 

The criterion for the equilibrium is therefore 

a^ A + bps + ■■■ = c^c + ■■■ , (87) 

which through Eq. Ij35|l leads to the relationship among 
the species abundances 



( Na\ d a 
\ NZ A 



NZ B 



NZ C 



We have developed a mathematical framework in terms 
of a Hclmholtz free-energy model which aims at de- 
scribing equilibrium properties of D-dimensional multi- 
component fluids. This theory is an extension of a former 
model devoted to one-component fluids and combines 
free-energy minimization methods with space partitions 
based on the available volume v to a particle and the oc- 
cupation number distributions of this variable. Thus we 
obtain an unified formalism that treats thermodynamics 
and structure of complex fluids simultaneously. 

The theory is presently developed for the case of di- 
lute fluids, where inter-particle interactions are evaluated 
from pair potentials. Two closure relations Ij23(l and lj38|l 
are required to obtain a closed formalism. The former is 
associated to the particle number conservation J2J and 
provides the probability that a particle had a specific 
species as the closest neighbor. The second closure re- 
lation l)38|l specifies the available-volume variable and, 
therefore, is linked to the filling space condition Q . Con- 
sequently, unconventional thermodynamic variables were 
derived by introducing space partitions in the thermo- 
statistical analysis of fluids. They are Lagrange's mul- 
tipliers related to the specific particle number and vol- 
ume conservations, which describe thus chemical poten- 
tials and the pressure, respectively. In addition, we have 
obtained general expressions for the excess chemical po- 
tentials and chemical equilibrium conditions of non-ideal 
fluids, as well as several pair distribution functions. 

The theory provides a straightforward way to calcu- 
late thermodynamic quantities and structure properties 
of gas mixtures. As illustrative examples, we have con- 
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sidered several hard sphere systems for which simple an- 
alytical expressions were derived. The equations of state 
obtained for pure hard-sphere systems are in good agree- 
ment with the scaled-particle and Percus-Yevick solu- 
tions valid for low-densities. We reckon that this agree- 
ment confirms a posteriori the validity of the several ap- 
proaches adopted. Besides, well-known solutions for van 
der Waals mixtures were recovered easily by using gener- 
alized Kac potentials. We also have derived closed ana- 
lytical forms for thermodynamic quantities of square- well 
mixture fluids. These evaluations constitute a reference 
for future studies of short ranged attractive systems. The 
influence of repulsive and attractive forces on the space 
partition among particles was analyzed. The results thus 
obtained express a close relationship between statistical 
geometry and statistical mechanics. 

The formalism introduced in the present work provide 
a useful tool to investigate the structure and the effec- 
tive interactions in chemical complex fluids. It is particu- 
larly well suited to study some astrophysical gases where 
detailed descriptions of species and interactions are re- 
quired. This occurs, in particular, in atmospheres and 
circumstellar envelopes, where the density can be hardly 
higher than some r\ = 0.01, so that the approximations 
for low densities given here can be used safely. We con- 
sider that present theory has other many applications in 
applied problems of current interest, such as studies of 
equilibrium properties of colloidal and protein solutions. 

Among the most important results obtained in this 
work, we can quote : 

• A functional (or series of functionals for mixing) 
contains the entire thermodynamical information of 
the system and leads to a classification of particles 
according to their distance to the closest neighbor 
through the assignment of a given volume available 
to each particle. So, it is thus possible to distin- 
guish the atoms according to the different degrees 
of perturbations carried by their neighbors. To our 
knowledge, there is no other model able to make 
this distinction up to now. 

• The current formalism is able to answer to the fol- 
lowing inquiries : 

— What is the probability of finding a particle of 
class j at a distance r from another particle of 
species i whose closest neighbor is of class k, 
knowing that r, i, j and k are arbitrary ? 

— Idem as before, but for a neighbor k situated 
at a given (arbitrary) distance from the parti- 
cle of class i ? 

• The chemical equilibrium of a fluid and the clas- 
sification of particles according to the distance to 
the respective closest neighbor, are done simulta- 
neously and self-consistently with the laws of ther- 
modynamics. 



• As far as we know, it is the first time that a theory 
is able to provide thermodynamical variables asso- 
ciated with the distribution of the space between 
the particles of different chemical species directly. 
The variables are related to the conventional ther- 
modynamical variables explicitly, such as pressure, 
chemical potential, etc. 

The description level attained by our formalism over- 
rides whatever known pre-existing model and constitutes 
an unparalleled tool to represent chemical mixing at 
low densities. It also provides conceptual elements on 
the connection between statistical mechanics and spatial 
statistics, which has never been presented up to now. 
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APPENDIX A: THERMODYNAMIC IDENTITIES 

In what follows we derive three thermodynamic iden- 
tities which enable us to derive relations l(T§|) . l(2T?|) and 
(|22|1 . The starting point is the equation that gives the to- 
tal particle population in terms of the occupation number 

N; 

in r • 



• Jo 



jmv dv . 



(Al) 



If in IjAljl we replace Nj mv by its expression given by 
(|14f> and derivate it with respect to a given population 
Nu by keeping constant the temperature, volume and the 
remaining populations (Nj m9 m), we shall obtain 



J2jm So Njmv 



_dj_ 



dN il 



dv 



v _dj_ 



(A2) 



dNu 



dv = 0. 



The second line in l|A2(l was obtained using relations © 
and 10}. We also note that 



J2jm Jo Njmv g 3 Nu dv 



and 



N, 



^J 3mv dN, 

jm 



jmv 2V 



3 dv = N, 



r^dv = e», 



dU vl 



3 ON, 



a 



(A3) 



-e l . (A4) 
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Eqs. ljA2j) to HA4|) lead to the identity 



V 



E^ 



dfjj 



dN a 



= . 



(A5) 



Then, Eq. Ijl9|) can be obtained easily from relations l|18(l 
and (|A5(1 . In a similar way, the partial derivative of Eq. 
(|A1|) with respect to the volume produces 



V f V AT I" ■ d7 dMjn, j de jmv de, 
Z^jmJo ^jrnv V gy gy T gy T gy 



V dV 



-Y N. 



jm dV 



dUi, 



dv 



dV 



with 



jm J 



de,- 



dV 



-dv. 



(A6) 



(A7) 



In Eq. I|A6(1 we have ignored the terms with Nj m ^ v=v ^, 
since they become zero at the thermodynamic limit. Ac- 
cording to Eq. lfT3|) , we have 



de,- 



E 



Nid- 



dV V 
In general dy^jmv.i — 0, so that 



~2V 



C — — ^ ^ / A ;„><•< iinrd'l' 



V 



(A8) 



(A9) 



Knowing that P = —dF/dV, relation l|20() is obtained 
straightforwardly from Eqs. l(T%|) and l|A10[l . 



The procedure followed to obtain relation (|A5|I can also 
be used to obtain the partial derivative of Eq. IjAlJ) with 
respect to the population of class i particles (iVj) by con- 
sidering constant the temperature, volume and the pop- 
ulations of the remaining species (Nj^i). We then use 



'. • E v ' ,Vv 



and 



de 



dNi 



2V 



The third identity that we can obtain is 



V 



ON, 



N 



-E 



j^i jm 



]m ON, 



ON,. 



(All) 



(A12) 



= 0. (A13) 



N4: 



Hence, from relation l|A6p we obtain the following iden- 
tity 



dv ^ J 



dV 



dV 



V 



= 0. (A10) 



Equation l|A13|) enable us to find an expression for the 
chemical potential ^ in terms of the thermodynamic po- 
tentials Hj m [Eq. P2*|]. 
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